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Abstract. We study relations between the eigenvectors of rational matrix functions on the 
Riemann sphere. Our main result is that for a subclass of functions that are products of two 
elementary blocks it is possible to represent these relations in a combinatorial-geometric way 
using a diagram of a cube. In this representation, vertices of the cube represent eigenvectors, 
edges are labeled by differences of locations of zeroes and poles of the determinant of our 
matrix function, and each face corresponds to a particular choice of a coordinate system on 
the space of such functions. Moreover, for each face this labeling encodes, in a neat and 
efficient way, a generating function for the expressions of the remaining four eigenvectors 
that label the opposing face of the cube in terms of the coordinates represented by the 
chosen face. The main motivation behind this work is that when our matrix is a Lax matrix 
of a discrete integrable system, such generating functions can be interpreted as Lagrangians 
of the system, and a choice of a particular face corresponds to a choice of the direction of 
the motion. 



1. Introductions 

Over the last 25 years a lot of research efforts have been directed towards the study of 
discrete analogues of integrable systems and, in particular, on how to adapt the existing 
methods and techniques from the classical theory of differential completely integrable sys- 
tems to the discrete case. The present paper is a small part of a larger project that aims 
to understand the Lagrangian structure of discrete integrable systems directly in terms of 
their Lax matrices and is motivated by work of Krichever and Phong [KP98], who obtained 
expressions for a universal symplectic form and elementary generating Hamiltonians on the 
space of Lax matrices for continuous completely integrable systems. Note that a univer- 
sal formula for a Lagrangian description of integrable systems in the continuous case was 
obtained by Zakharov and Mikhailov [ZM80]. 

For discrete completely integrable systems, the Lagrangian point of view has been de- 
veloped in the classical papers by Veselov [Ves88, Ves91] and Moser and Veselov [MV91]. 
Recently a very promising approach to the study of Lagrangian structure of integrable lat- 
tice equations in terms of Lagrangian multiforms has been proposed by Nijhoff and Lobb 
[LN09]. This approach is related to the notion of multidimensional consistency formulated 
by Bobenko and Suris [BS02] and independently by Nijhoff [Nij02]. Although our approach 
is much more elementary and is based on the notion of the refactorization transformations, 
as in the original work of Veselov and Moser, the resulting combinatorial diagrams encoding 
such transformations are reminiscent of the multidimensional consistency approach, and it 
will be very interesting to see if this is more than a pure coincidence. 
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The appearance of refactorization transformations in the theory of discrete integrable 
systems is not very surprising. Indeed, a discrete analogue of the Lax Pair representation is 
the isospectral transformation L(z) = R(z)L(,2)R(,s) _1 , where L(^) is the Lax matrix of the 
system and R(z) is the evolution matrix that has to be chosen in a special way dependent on 
L(z). When ~L{z) is a rational matrix function with the fixed singularity structure, one can 
specify elementary evolution matrices R(.z) using pairs of points p± on the spectral curve T of 
L(z) and the eigenvectors of L(z) at those points, and choosing those points to lie above the 
points of the determinantal divisor of detL(z) corresponds to factoring L(z) = L 1 (z)B 2 (-z) 
and mapping it to L(z) = B 2 (-2)Li(;z), i.e., by choosing R(z) = B 2 (^). 

In describing such refactorization transformations it is essential to choose a good coor- 
dinate system on the space of Lax matrices. The natural candidates for the coordinates 
are the eigenvectors of L(Y) and L(z) _1 . However, for the refactorization transformations 
another important class of vectors consists of vectors defining elementary building blocks. 
The relationship between these two sets of vectors is quite complicated, and in attempt 
to understand it we noticed a very elegant way to represent it using geometric diagrams 
described in the present paper. Although we restrict our attention to the quadratic case 
when L(z) is a product of two blocks and when there is no essential difference between the 
vectors of the additive and multiplicative representation, it is still quite interesting. Indeed, 
such refactorization transformations describe, for example, the change of polarization in the 
interaction of soliton solutions of the matrix KdV equations [GV04] and integrable discrete 
vector nonlinear Schrodinger equations [PAB06]. Further, in [Dzh09] we showed that for 
two-dimensional matrix functions the directional dynamics of the eigenvectors is described 
by difference Painleve-V equation. 

We plan to use this tools developed in this paper to study the higher dimensional case of 
products of three and more elementary blocks and, in particular, the relation to Yang-Baxter 
maps, in a separate publication. 

This paper is organized as follows. In the next section we briefly describe a particular 
class of Lax matrices that we consider, their additive and multiplicative representations and 
the notion of an elementary divisor — a building block for a multiplicative representation. 
In Section 3 we explain a visual representation of some linear equations involving elementary 
divisors and in Section 4 we use this visual representation to describe the relations between 
eigenvectors of a quadratic Lax matrix, and how to use the resulting cube diagram to generate 
the Lagrangian functions of the refactorization transformations. 



The main goal of this section is to describe the space of Lax matrices, their additive 
and multiplicative representations and related coordinate systems, and briefly review the 
relationship between isomonodromic transformations and discrete Painleve equations, see 
[Dzh09] for details. 

2.1. The space of Lax Matrices. We consider the space C of rational mxm matrix func- 
tions L(z) satisfying the following conditions. We assume that L(z) is regular, diagonalizable 
(and diagonalized) at some normalization point zq that we take to be z® = oo, 
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(2.1) 
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We also assume that the singularity structure of L(z), i.e., points where L(z) has a pole or 
becomes degenerate, is accurately represented by its determinant, and that the determinant 
is generic, i.e., it has only simple zeroes and simple poles. This amounts to requiring that 
L(z) is holomorphic except for simple poles at the points zi,...,Zk, M(z) = L(z)- 1 is 
holomorphic except for simple poles at the points Ci> ■ ■ ■ > Cfc> an z i an d Cj are distinct, and 

II a (* - Ca) 



(2.2) 



det L(z) = pi • • • p % 



Uk( z - z k) ' 

These conditions mean that the residues L, = res Zi L(z) and Mj = — res^. M(z), where 
the negative sign here is just for convenience, are matrices of rank one. Using the f symbol 
to indicate a row vector, we have: 

k 



(2.3) 
(2.4) 

(2.5) L(z 



Uz) 



detL(z) 



1=1 



Z - Zi 



where L r 



diag{p 



i- 



,p m } and Li = a^b], 



Pi ■ ■ -Pm^k 



i=i\ z z h 



-i 



M{z) 



where M r 



T - 1 



M, 



c d f 



The above representations of L(z) and M(z) are called additive representations and the 
vectors aj, b\ (resp. Cj, dj) additive eigenvectors of L(z) (resp. M(z)). Note that these 
eigenvectors are also characterized by L(£ a )c a = M(2j)aj = and d£,L(C a ) = b}Mj = 0. 

Let V = J2i z i ~ J2i d be the divisor of the (determinant) of L(z) and denote the space of 
matrices L(z) satisfying conditions (2.3)-(2.5) by Mf. Then eigenvectors give coordinates 
on this space, as described by the following lemma. 

Lemma 2.1. Generically, the collection {a^, dj c }^ =1 (or the collection {cfc,b^,}^ =ly ) is a co- 
ordinate system on the space Ai® . 

Proof. Consider the equations M(zfc)a^ = and djL(£j) = 0: 



(2.6) 

Then if the matrix 
(2.7) 



L 



" d! 

c 



&k 



z k — Ci 

is invert ible, 



0, djL + 



k=l 



_dja^ 1 

Ci ~ z k 



0. 



L 1 a^. 



z k — Ci 



z k — Ci 



dL . 



□ 



2.2. Elementary Divisors and Multiplicative Representations. To define a multi- 
plicative representation of L(z) we first define its building blocks. These are rational matrix 
functions of the special form 

Q 

(2.8) B(z) = I H , where G = fg^ is a matrix of rank 1. 

Z — Zq 



We call such matrices elementary divisors [Dzh09]. If tr(G) = g^f 7^ 0, we can consider 
instead of G a rank-one projector P = f(g^i)~ l g^ . Then we can write 



(2.9)B W =B( Z ;,„,C„)=I + f^, d*B W = i^, B W - = I+f^M. 

Z — Zo g'l z — Z z — (,0 g 1 1 

If we now pair the zeroes and poles of det ~L(z) in some way as (£ s , z s ), then, for any such 
pair, there is a multiplicative component of L(z) of the above form. There are two ways to 
think about such components — we can look at factors or at divisors, see also [Bor04]. 

Definition 2.2. We say that elementary divisors B^z) (resp. B l s (z)) corresponding to pairs 
(Cs,z s ) are right (resp. left) divisors of L(z) if L(z) = L£(z)B'£(z) (resp. L(z) = B z s (,z)L/ s (z)) 
where L^z) (resp. L^(^)) is regular at z s . Further, we say that elementary divisors B s (z) 
corresponding to pairs (( s , z s ) are the factors of L(z) if h(z) = L Bx(z) • ■ ■ B k (z). 

In [Dzh09] we showed that left and right divisors can be written explicitly in terms of the 
eigenvectors of L(z), 

(2,0) B; W =i + ^g, Bl W = I + ^g, 

and so the coordinate systems described in Lemma 2.1 are just parameterizations of Lax 
matrices by left (resp. right) divisors. When L(z) has only two factors, which is the quadratic 
case that we focus on in this paper, there is no essential difference between divisors and 
factors. 

2.3. Re-factorization transformations and d-Py. We now consider the refactorization 
transformation in the quadratic case. Let 

(2.11) L(z) = L B 1 (^)B 2 (z) = B< (z)L B^) = B[(z)L B r 2 (z) 
Then we can consider either isospectral transformation with R(z) = B\(z), 

(2.12) L(z) = B l 2 (z)L B{(z)^ L(z) = BJ(z)Bj,(z)L = B^)L B^), 
or isomonodromic transformations 

(2.13) L(z) = B l 2 (z)L B[(z)^ L(z) = B[(z + l)B l 2 (z)L = B l 2 (z)L B[(z) . 

When m = 2, the resulting phase space is two-dimensional and it possible to introduce the 
so-called spectral coordinates (7, 71") on the space of Lax matrices. Then, in the isomonodromic 
case, the spectral coordinates of L(z) and L(z) are related by the difference Painleve V 
equation in the Sakai's classification [SakOl], see [Dzh09]: 

(2.14) 7 + 7 = » + (2 4- Mh ~ 3 + b) + Mh ~ Zl + C ' + 1} . 

It- pi 7T - p 2 

(2-15) 7T7T = PlP2 - — j-. 

(7 - zi)(7 - Ci) 
In [Dzh08] we proved the following 



Theorem 2.3. The equations of both the isospectral and isomonodromic dynamic can be 
written in the Lagrangian form with 

(2.16) £(X, Y,t) = (z 2 - z x {t)) log(x{x 2 ) + ( Zl (t) - log^I^x,) 

(2.17) + (Ca - Ci(0) log(ylLo 1 y 2 ) + (Ci(*) - z*) log(x{y 2 ), 

where X = (x^xl,) andY = (yi,y|), in the isomonodromic case z\{t) = z\—t, Ci(£) — Ci — t, 
and in the isospectral case Z\(t) = Z\, (i{t) = Ci and C(X,Y) is time-independent. 

Here xj, Xj, yj, y 3 - are some eigenvectors of L(z) and L(z). 

The main goal of the present paper is to give a combinatorial-geometric representation of 
this generating function. 



3. Linear Equations given by Elementary Divisors and their Geometric 

Representation 

In studying the relationship between various eigenvectors of L(z) we mainly work with 
linear equations of the form v ~ B(z*)w and ~ w*~B(z*), where ~ means that the two 
vectors are proportional and H(z*; z ,(o) is a matrix of the form (2.8) evaluated at some 
point z* . In this section we give a diagrammatic representation of such equations using the 
language of elementary triples. 

3.1. Basic Definitions. In what follows, let V be an n-dimensional complex vector space 
whose elements are column vectors, V ~ C n , and let V' be its dual-space. We think of 
elements of as row-vectors and denote its elements as q} G ~ (C n y. Our main 
definition is the following. 

Definition 3.1. Let \ G PV, % = 1, ... 3, be one-dimensional linear subspaces of V and 
//' G PVt be a one-dimensional linear subspace of VL We say that (Ai,A2,As) form a 
-based triple with parameters (ai,a 2 ,a 3 ), a, G C\{0}, if the following conditions hold: 

• Ai ^ ker ^ for all i; 

• Aj 7^ \j for % 7^ j but Aj C Xj + A& for all i,j, k. 

3 

• Qi A = 011 A | 012 \ | ft3 x =0 

1 A* f Ai 1 /utA 2 2 /itA 3 3 

The last equation above adjusts the "slope" of \ w.r.t. Xj and A& and it should be inter- 
preted in terms of spanning vectors. To that end, let Aj = Span c {pj} and jj) = Span c {q"l"}. 
Then we require that the vectors pi, p 2 , and p 3 are linearly dependent but pairwise inde- 
pendent, qtpj 7^ for all i, and the last equation takes the form 

3 

\ QJj Oil a 2 a 3 

(3-1) > Pi = "j - P 1 + ~^P2 + ~^P3 = 0. 

^qt Pi qt pi qtp 2 q tp 3 

Note that this equation is homogeneous and so is independent of the choice of basis vectors. 
We represent such triples by diagrams of the form 

5 



P3 




where we use the black circles for subspaces of V (or their basis vectors) and white circles for 
subspaces of V' . In what follows, we switch freely between subspace and vector formulations. 

It is clear that a^s are defined up to a common multiplicative constant and satisfy 

(3.2) ai + a 2 + a 3 = 0. 

Moreover, given qt, there is a unique vector [pj] = [pi] q t G A, normalized by the condition 
q f [p] = !> i- e -> 

(3.3) [ Pi ] = ~r~ ^ A*, 

q T Pi 

where p» is any non-zero vector in Aj. Then, in terms of the normalized vectors, the linear 
dependence equation (3.1) becomes 

(3.4) ai[pi] + a 2 [p 2 ] + a 3 [P3] =0 or [p;] = a ^ P ^ +Qfc [P fc ] ; ^ y. & \\ distinct, 

ciij + a k 

which explains the interpretation of this equation in terms of "slopes" . 

Definition 3.2. We define the X-based dual triples (fit, fi 2 , ^3) with parameters (/3i, (3 2 , /3 3 ), 
Pi G C\{0}, where A G PV is a one-dimensional linear subspace of V and //J G PV^, 
% = 1, ... 3, are one-dimensional linear subspaces of in exactly the same way, and represent 
them by the following diagrams: 



(4 q 3 




Finally, the behavior of the triples under the action of GL(V) is described by the following 
obvious Lemma. 

Lemma 3.3. Let A G GL(V). 

(i) If (Ai, A2, A3) is a -based triple with parameters (oti, a. 2 , ot$) , then (AAi, AA 2 , AA3) is 
a /J A~ 1 -based triple with the same parameters (ai, a 2 , 0(3); 

(ii) dually, if(fi\, /4> /4) be a X-based triple with parameters (flx, f3 2 , (3 3 ), then (/i{A _1 , /ijjA -1 , /igA -1 ) 
is an AX-based triple with the same parameters (0%, {3 2 , /? 3 ). 

6 



3.2. Elementary Divisors and Triples. We now turn our attention to the study of the 
relationship between triples and elementary divisors. The following Lemma is immediate. 

Lemma 3.4. Let 

3.5 B = I+- =1+ pq T =I+— pq T , 

z* — z qTp z* — z z* — z 

[p] = [p]qt an d [q'] = [q ]p- Note that we can normalize either p or q' ; but not both. The 
choice of which one should be normalized depends on whether we think o/B as acting on the 
elements of V or . Then 

(i) Vectors p and qt are the eigenvectors o/B, 

(3.6) q^B = — q*, Bp = — p, where detB = — . 

Z — Zq Z — Zq Z — Zq 

(ii) For any vector w such that q^w ^ 0, and for any row-vector such that w^p ^ 0, 
using normalized vectors we get 

(3.7) (B - I)[w] qt = ^[ P ] q t, [wt] p (B - I) = ^^[qt] p . 

Z — Zq Z — Zq 

(Hi) Ifv = Bw, then [v] q t = (z* — zo)/(z* — Co)B[w] q t, and so (w,p,v) /orm aq^-triple with 
parameters (z* — z , z — ( , ( — z*). Similarly, ifx^ = w^B, then [v^] p = (-2* — zq)/(z* — 
Co)[ w ^]pB, and (w*, p^, v*) form a p-triple with parameters (z* — z , z — ( , ( — z*). 

(iv) Conversely, for any q} -triple (pi,P2,P3) with parameters (a±, a 2 , ol?), 

(3.8) [p 3 ] = (l + -[p 2 W) [pi] = —J— /x (ailog(qt pi ) + a 2 log(qtp 2 )) , 

«i + «2 \ ai J 011 + a 2 Cq T 

and for any p-triple (q{,q2,q3) with parameters (fa, fa, fa), 

<V» M = ATft [q!l ( r + I p|qSl ) = ATA | ( A log(qIp » + * log(qSp) ) ' 

3.3. Gluing Properties. The following two Lemmas show that any two triples with match- 
ing parts can be "glued" together and then completed into a cube in such a way that any 
vertex of the cube makes a triple. More precisely, the Face Lemma shows how to glue two 
triples of the same type that share two vertices, and the Edge Lemma shows how to glue 
two triples of different type if the center of one is the end of the other and vise-versa, and 
the common edge parameters coincide. 

Lemma 3.5 (The Face Lemma). Let the pairs pi,P2 G V and q\,q\ G be linearly 
independent and not-orthogonal, qjpj 7^ 0, and let ot\, a 2 , fa, fa be any non-zero complex 
parameters. Then the quad 



can be completed to a cube with 

(3.10) a 3 = -(a 1 + a 2 ), 03 = ~(0i + #2), 7i = -(<*i + A), 72 = -(a 2 + /3 2 ), 
and some oi/ier parameters a 4 , a 5 , /3 4 , /3 5 , 




pi 



71 





a 2 




ft . 








.,-'"4 






1 


05 



4 



P2 



ft 



72 



Q 2 



ft 



50 that at each vertex the corresponding triple condition is satisfied. Moreover, if the param- 
eters satisfy the additional cyclic condition 

(3.11) ai + a 2 + ± + 2 = 0, 

on the initial face, then the same condition is satisfied on the remaining five faces, and so 

(3.12) a 4 = —0 2 , «5 = ~0i, 04 = -«2, 05 = -Oii- 
Further, let 

(3.13) C = ai log(qjpi) + a 2 log(qjp 2 ) + X log(q 2 pi) + 2 log(q t 1 p 2 ). 

Then 



(3.14) [P!] q t=«3g, 



[Ql] P i = 7i 



<9£ 
9pi 



[Q 2 ] P2 = 72 



9pi 



Proof. First, note that the normalized outside vertices are uniquely determined by the given 
data, 



(3.15) 
(3.16) 



«l[Pl] q t +«2[p 2 ] q t 
(Xl + «2 

rn ti _ gjJgiU+Abglpi 

CKi + pi 



A[pi] q t +/3 2 [ P2 ] q t 



/3i + /3 2 
rn ti _ Q 2[qI] P2 +/3 2 [q 2 ] P2 

[^2JP2 — I 

a 2 + p 2 



Thus, to prove the theorem we need to show that it is possible to find the parameters 

«4, «5, /?4, 05 so that 

(3.17) a A + a 5 = a x + a 2 , 0a + 05 = 01 + 02, «4 + 0a = "1 + 0i, "5 + 05 = "2 + 02, 
and that the following equations hold at the outside nodes: 



(3.18) 
(3.19) 
(3.20) 
(3.21) 



qS 

Pi 
p 2 



«4[Pl] Q t+/3 4 [P2] Q t +7i[Pi] Q t =0 

«s[Pi] Q t +/3 5 [P 2 ] Q t +7 2 [p 2 ] Q t =0 

«4[QI]pi + «5[Q 2 ]Pi + «3[ql]Pi = 

/34[Ql]p 2 +/35[QS]p 2 +/3 3 [qS]p a = 0. 



Using the change of normalization rule 

(3-22) [p] qi = ^[P] q t, 

the first equation, when decomposed in the spanning basis pi,P2, becomes the following 
system: 

"4 qlPi \ Qlpi ^ ( & 4^2 \ R Qlpi l r 

Oil— e h ^— 5- — f— • Pi — e — = ai + pi 



(3.23) 



tti + ^QlPi/ q f lPl \/3i + I3 2 Q\P 2 J qt pi 

«4 qlPl \ Q2P2 . ( P4 q^2 , , 



«i + «2 QjP! y qip 2 v & + & q1p 2 



Solving it gives 

(3.24) a, = /3 2 (a 1 + a 2 )(a 1 + /3 1 )[Ql] Pl [P 1 ] qt ^Pi^Pi^ 

41 a 1 /3 2 q 2 PiqiP2 - a 2 /3iq[piq 2 p 2 

(3.25) /3 4 = a2 (0 1 +0 2 )( ai +A)[Q t 1 ] P1 [P2] q t - -qjp^qjp 2 , 

2 «i/?2q 2 PiqiP2 - « 2 /3iq|piq 2 p 2 

which, together with the expressions 

,a 9fi ^ r n ti rp 1 _ zl ~ ^P^^ ~ t^/WiPiq^ 

(3.2b) NdpdP^t - ^ + Q2)(ai + ft) • qtp iq tp iq tp 2 

/o ova rnti \v> 1 - 1 «i/32q^Piq t iP 2 - (gg ~ ftQ/WiPiq^ 

W-^'J Nl Pi *2 Lt - / fl . fl \ / ,~TT ' + 1 1 ' 

q2 (Pi + p 2 )(ai + Pi) qlpiq 2 piq 2 p 2 

where k = a\ + a 2 + p\ + /3 2 , gives 

«i/32q 2 PiqIp2 



(3.28) a 4 = — /3 2 ~\~ k 



ai/3 2 q 2 Piq t iP2 - Ofe/WPiQ^ ' 
-a 2 /3iq t 1 piqjp 2 



(3.29) h -n, • /,• . 

ai/32q 2 PiqiP2 - a2piq|piq 2 P2 

and so it is easy to see that 014 + p4 = k — o; 2 — p2 = cti + p\, as required. Similarly, using 

/o qnA r n ti rp 1 1 ai^gPiqiPg - a 2 (/3i - ^qlpiggPg 

(d.duj [^W^iLt - 7- — —\r —q—\ ' ; 1 1 

q i (ai + a 2 )(a 2 + /^ 2 ) qipiqip 2 q 2 p 2 

/o qi\ rnti [pi ~1 (^x - /^qjpiqlpg - c^/WiPiq^ 

^d.di; [ t 4 2 J P2 [^2j t - 775— tw- — ~o\ + t t ' 

Q2 (Pi + P2)(«2 + p 2 ) q 2 piqip 2 q 2 p 2 

and the equations at the node Q 2 gives 

( oo \ n a , „ -« 2 Aq t iPiq 2 p 2 

(3.32) a 5 = —pi + k 



(3.33) (3 5 = —ai + k 



ai/3 2 q 2 piq t iP 2 - a 2 /3iq t iPiq 2 p 2 ' 
«i/3 2 q 2 piqlp 2 

ai/3 2 q 2 Piq t iP 2 - « 2 /3iq t iPiq 2 p 2 

9 



Similar computations show that equations at the nodes Pi and P2 give the same values of 
the parameters, and the other constraints are easy to check. □ 

Lemma 3.6 (The Edge Lemma). Let (pi, P2, P3) be a q^-triple with parameters (ax, a 2 , 0:3), 
(Qi) ^3) a P3-triple with parameters (Pi, (3 2 , P3) with parameters /3 2 , P3), and let 
"3 = h, 



93 



4 



<*3 = fa 




Then it is possible to glue the triples along the p^^-edge and, by assigning parameters 71 
and 72 along the piqj and p 2 <^ 2 edges so that the sum of parameters in the new quads is zero 
(i.e., 71 = «2 — Pi and 72 = a± — $2), complete the resulting "butterfly" configuration to a 
consistent cube, 




Proof. After choosing 71 and 72 according to the above rule, this follows immediately from 
the Face Lemma 3.5. 

□ 



4. Eigenvectors of Quadratic Lax Matrices and Refactorization 

In this section we show that a cube diagram described in the previous section can be used 
to represent the relationship between eigenvectors of quadratic Lax matrices. Also, we show 
that the three axes of symmetry of the cube connecting the centers of opposing sides can be 
give the following interpretation. One of the axis corresponds to switching the coordinate 
system on the space from left to right divisors. The other two axes correspond to two 
different directions of the refactorization dynamics corresponding to two possible pairings 
between two zeroes and two poles of det L(z). For each of this transformations we produce 
the explicit formula for the generating function using the labeling of the cube. 

4.1. Relations Between Eigenvectors of Quadratic Lax Matrices. Consider a qua- 
dratic Lax matrix ~L(z) given in the additive form as 

(4.1) L(z) = L + ^- + ^L, L(z)- 1 = L 1 C « d « - C ^ 



(4.2) detL(z) = detL. 



Z - Z 1 Z- Z 2 ' Z - (a Z-C/3' 

(Z ~ ( a )(z - (p) 



zi)(z - Z 2 ) 
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Pairing z\ with ( a = d and z 2 with = ( 2 we can write L(z) and L l (z) in the multi- 



plicative form using factors Bi(z) = I + 
(4.3) L(z)=L B 1 (z)B 2 (z), 
and so Bi(*) = L^B^L, 



z-c,i p»q 

z—z. 



L^)- 1 = B 2 (^)- 1 B 1 (^)- 1 L- 1 , 



I + Zl - (a L ° la t ld " L ° and B 2 (z) = B r Jz) = I + This 

z-zi dLai Zy 1 2W z-z 2 blcp 

can also be seen directly by looking at the residues of L(z) at zi and 22 and of L(z) _1 at Co 
and (p; in fact, that is how the expressions for the right and left divisors are obtained. The 
same residues also give the following collection of equations: 

(4.4) dlh B 2 ( Zl ) ~ b\, B!(z 2 ) C/3 ~ L B 2 (Ca)c a ~ L^ai, d^B^) ~ b 2 . 
Representing these equations using elementary triples and gluing, we get the following result. 
Theorem 4.1. The eigenvectors ofL(z) are related by the following cube diagram: 



d^Lo 



22 - da 



L 1 



a i 





z l — Cq 


<^ 

1 


CM 

M 
1 

CH- 


1 


Ca ~ z l 





J a 2 



I 



')9 



Note that in this diagram, the top (resp. bottom) faces of the cube correspond to parame- 
terizing the space of quadratic Lax matrices by left (resp. right) divisors, and so the vertical 
axis of the cube correspond to the change of coordinates between these two different systems. 
Two horizontal axes correspond to the refactorization dynamics. 



4.2. Refactorization Dynamics. Consider now an isospectral dynamic L(z) h-> L(z) 



H,(z)'L(z)'R(z)~ 1 , the isomonodromic case is similar. If we take R(z 
where the notation B 

this becomes a refactorization transformation 



B 2 (z) = WUz) 



,8;2V 

p. 2 \z) explicitly specifies the zero and the pole of the elementary divisor, 



(4.5) 



L(z) = L B a;l (z)Bp. 2 (z) h-> L(z) = B /3;2 (z)L B Q;1 (z) = L B a;1 (z)B /3 . 2 (^) 



Since we can think of the refactorization transformation as switching the roles of the left 
and right divisors, B^ 2 (z) = B r ^ 2 (z) = B l ^ 2 (z) and B Qjl (z) = L^B^OJLq = B£ (1 (s), 
using (2.10) we get the identifications a 2 = eg, d^ = b 2 , c a = L 1 a 1 , and b{ = d^L . Thus, 
on the diagram below labels on the back face correspond to the coordinates of L(z), labels 
on the front face correspond (after twisting a 2 by Lq 1 and dt by L ) to the coordinates of 
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and the generating function encoded by the right face of the cube is, up to some minor 
change of notation, exactly the Lagrangian from Theorem (2.3), 

(4.6) £((a 2 , b{), (a 2 , b{)) = (z 2 - z x ) \og{h\k 2 ) + (z x - Cp) log(b}a 2 ) 

+ (C/3 - Ca) log(blL x a 2 ) + ((a ~ z 2 ) log(bjLo 

Note that using the generating function given by the left face of the cube corresponds to the 
backwards motion generated by R(z) = B^. x (z). 

The remaining axis of the cube corresponds to pairing z\ with and z 2 with £ w . In other 
words, taking R(z) = B£. 2 (z) (or R(z) = BL(z) for the backwards motion), 

(4.7) L(z) = L B /3;1 (z)B a;2 (2;) ^ L(z) = B a . 2 (z)L B P;1 (z) = L B ftl (^)B a . 2 (^, 
we get the following map 




and its Lagrangian generating function 

(4.8) £((a 2 , b|), (a 2 , bj)) = (z x - z 2 ) \og(h%) + (C« - z x ) log(b t 1 a 2 ) 

+ (^ - Ca) log(blL x a 2 ) + (z 2 - Cp) log(bjL- x a 2 ). 

5. Conclusions 

We explained a neat and efficient way to encode the structure of refactorization transfor- 
mations for quadratic Lax matrices and their generating functions using cube diagrams. This 
approach also gives different ways of choosing coordinate systems on the space of such matri- 
ces — each choice corresponds to a face of the cube. It would be very interesting to see if this 
approach can be generalized to Lax matrices with more than two factors. Since, in view of 
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the gluing properties, cube diagrams are rigid, we expect it to results in higher dimensional 
configurations relating such cubes, where each cube represents a particular transposition of 
factors. In particular, the case of three factors is related to the structure of Yang-Baxter 
maps and we plan to consider it in a separate publication. 
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